arXiv: 1501.02817v2 [cond-mat.str-el] 26 Apr 2015 


Multiscale Entanglement Renormalisation Ansatz for Spin Chains with Continuonsly 

Varying Criticality 


Jacob C. Bridgeman, Aroon O’Brien, Stephen D. Bartlett, and Andrew C. Doherty 
Centre for Engineered Quantum Systems, School of Physics, 

The University of Sydney, Sydney NSW 2006, Australia 

We use the multiscale entanglement renormalisation ansatz (MERA) to numerically investigate 
three critical quantum spin chains with Z2 x Z2 on-site symmetry: a staggered XXZ model, a 
transverse field cluster model, and the quantum Ashkin-Teller model. All three models possess a 
continuous one-parameter family of critical points. Along this critical line, the thermodynamic limit 
of these models is expected to be described by classes of c = 1 conformal field theories (CFTs) of 
two possible types: the S^ free boson and its Z2-orbifold. Our numerics using MERA with explicitly 
enforced Z2 x Z2 symmetry allow us to extract conformal data for each model, with strong evidence 
supporting the identification of the staggered XXZ model and critical transverse field cluster model 
with the S^ boson GET, and the Ashkin-Teller model with the .Z2-orbifold boson GET. Our hrst two 
models describe the phase transitions between symmetry protected topologically ordered phases and 
trivial phases, which lie outside the usual Landau-Ginsburg-Wilson paradigm of symmetry breaking. 
Our results show that a range of critical theories can arise at the boundary of a single symmetry 
protected phase. 


I. INTRODUCTION 

Quantum spin models are an area of extensive theoret¬ 
ical and numerical study, due to their relative simplicity 
and wide descriptive power. Simple spin models can ex¬ 
hibit a variety of exotic ground states properties, such as 
topological order^ and symmetry-protected topological 
order^. Studying the critical behaviour that describes 
a system at the transition point between two quantum 
phases can also be investigated with quantum spin mod¬ 
els but brings new challenges. At these critical points, 
the models become gapless and many of the exact results 
that have been proven for gapped systems break down. 
In particular, for critical systems the correlation length 
in the ground state diverges and the area law for en¬ 
tanglement entropy is violated^. In the thermodynamic 
limit, the behaviour at the critical point is described by 
a conformal field theory (CFT)^’®. 

A wide range of numerical methods have been devel¬ 
oped to study the behaviour of quantum spin chains, in¬ 
cluding their critical behaviour, despite the difficulties in 
working with a Hilbert space that grows exponentially 
in the chain length. Numerical methods based on tensor 
networks®’^ have been highly successful in recent years as 
an efficient method for studying a wide range of spin lat¬ 
tice models. In particular, for one dimensional gapped 
systems, methods based on matrix product states are 
known to efficiently describe ground state properties®. 
However, this behaviour does not generalise to gapless 
systems® (although c.f. Ref. 10). For a tensor network 
description to naturally describe a critical spin chain, 
is should capture the area law violation, and one such 
description is the Multiscale Entanglement Renormalisa¬ 
tion Ansatz^^ (MERA). The MERA has previously been 
shown to accurately and efficiently reproduce conformal 
data of critical spin chains^®^^^. 

In this paper, we use the MERA to numerically study 


three related spin models along a line of criticality: a 
staggered XXZ model, a transverse field cluster modeb®, 
and the quantum Ashkin-Teller modeb®d7 These mod¬ 
els are of interest because they have critical exponents 
that vary continuously as a function of model parame¬ 
ters. All three models possess an on-site Z 2 x Z 2 sym¬ 
metry that will play an important role in our analysis. 

For the first two models of interest, the critical line 
corresponds to a phase transition between a phase pos¬ 
sessing nontrivial symmetry protected topological or¬ 
der (SPTO) for this symmetry group^® and a trivial 
phase. Such systems are also of particular interest 
due to their connection to measurement-based quantum 
computation^®’^®’®®. There has been extensive investi¬ 
gation into symmetry protected topological phases, but 
only a few studies have investigated transitions out of 
them®^’®®. Such phase transitions are not part of the 
usual Landau-Ginsburg-Wilson paradigm of phase tran¬ 
sitions because there is no broken symmetry on either 
side of the transition. Despite the fact that the mod¬ 
els we consider here involve a transition between a fixed 
SPTO phase and the trivial phase, the boundary is de¬ 
scribed by a range of critical theories. 

In contrast, the Ashkin-Teller model possesses a con¬ 
ventional symmetry broken phase and for that model the 
phase transition of interest is between the symmetry bro¬ 
ken and trivial phases. Again, a range of critical theories 
describe the phase boundary. 

The thermodynamic limit of these three models at 
their critical points are described by a special class of 
CFT, namely those with central charge c = 1. These 
CFTs do not correspond to any member of the series of 
so-called minimal models that have been the main fo¬ 
cus of numerical research that has used MERA to study 
CETs^®’^^. This class of GET is characterised by a renor¬ 
malisation group with an exactly marginal field, leading 
to a single parameter family of theories that reproduces 
the continuously varying critical exponents of the under- 
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lying spin models. By varying the coupling in the spin 
models, we can tune the system along this line of criti¬ 
cality. This continuous variation of scaling dimensions is 
a distinctive feature of c = 1 CFTs^^. 

Here, we show that the MERA allows us to extract 
conformal data from all three models described above 
possessing a continuous line of criticality, and show that 
these data agree with the expected c = 1 CFTs includ¬ 
ing the variation of the conformal data along the line of 
criticality. Specifically, the staggered XXZ model and 
the transverse field cluster model are both shown to be 
consistent with the c = 1 5^ free boson CFT, whereas 
the quantum Ashkin-Teller model gives data consistent 
with the c = 1 orbifold CFT. Our results provide further 
evidence of the usefulness of MERA in describing critical 
systems in one dimension. 

This paper is structured as follows. In Section II, we 
review CFTs as they pertain to the three spin models 
studied. In Section III we review the details of the MERA 
algorithm. In Section IV we present results from the 
numeric simulations for each model, and compare with 
the expected CFT data. 


II. c = 1 BOSON CFTS AND ASSOCIATED SPIN 
MODELS 

We will consider three different one dimensional quan¬ 
tum spin models at their respective critical points. In 
each of these models there is a line of critical points along 
which the critical exponents vary continuously. These 
critical points are believed to be described by certain con¬ 
formal field theories with marginal operators that gener¬ 
ate flow along these lines of criticality. In this section, 
we briefly review these conformal field theories and in¬ 
troduce our three spin chain models. 


A. c = 1 CFTs 

Conformal field theories (CFTs) are field theories 
that are invariant under all conformal (angle preserv¬ 
ing) space-time transformations. They describe the ther¬ 
modynamic limit of critical lattice models^’®. Due to 
the abundant symmetry, a (l-l-l)D CFT is completely 
specified by knowledge of (*) its central charge c; {ii) 
the primary fields cj) and their scaling dimensions (h, h), 
which are eigenoperators of the scaling transformation 
with eigenvalues A — h + h; and (Hi) the operator prod¬ 
uct expansion coefficients for these fields. It is the pri¬ 
mary field scaling dimensions A that predict the critical 
exponents of the associated critical lattice model. 

We focus our attention on CFTs with c = 1, which 
mainly fall into two categories: the free boson theory and 
its Z 2 -orbifold^’®. These c = 1 boson CFTs are charac¬ 
terised by the presence of an exactly marginal primary 
field, with scaling dimension 2, meaning it is fixed un¬ 
der rescaling. Perturbing a CFT by an exactly marginal 



Figure 1 . Illustrative phase space diagrams and renormali¬ 
sation group (RG) flows for models with a unique RG fixed 
point (left) and a line of criticality (right). RG relevant opera¬ 
tors have A < 2 and grow under renormalization, whilst those 
with A > 2 are called irrelevant, since a CFT deformed by 
such an operator flow back to the fixed point. Deformations 
by an exactly marginal operator lead to a new conformal fixed 
point. As the size of the deformation is varied, the conformal 
dimensions vary continuously. 


term gives another CFT, leading to continuously varying 
families of theories (see Figure 1). This phenomenon is 
not seen in unitary CFTs with c < 1. These two continu¬ 
ous families of CFTs will correspond to the critical lines 
in our various quantum spin chains. 

The first relevant family of CFTs, the 5”^ boson CFT, 
is the theory of a massless free boson (p{x). We will be 
interested in periodic boundary conditions and we will 
choose a parameterisation such that x G [0,27r). We 
will focus on a compactified version of the free boson 
where (j) itself takes values on a circle of radius Rc, so 
that ip{x) = ^{x) -I- 2t:R c- As a result when we move 
around a circle in the x co-ordinate, Lp does not need to 
be strictly periodic but may twist around m times so that 
ip{x + 2 tt) = Lp{x) -\- 2TrmRc- The primary held scaling 
dimensions are known for this theory and generally have 
a non-trivial dependence on Rc- The free boson has a 
natural internal SO(2) symmetry given by the transla¬ 
tion (p{x) = ip{x) + 6Rc, where 9 G [0, 27r), as well as a 
Z 2 symmetry given by p>(x) —>■ —ip(x). 

The orbifold boson CFT once again involves a bosonic 
held (p{x) compactihed on a circle of radius Rq, so that 
(p{x) = (p{x)+2TrRo^’^ ■ The theory has two sectors. The 
hrst corresponds to the subspace of symmetric states of 
the free boson that map to themselves under Lp{x) —>■ 
—(p{x). The second sector corresponds to quantising the 
free boson with “twisted” boundary conditions for which 
(p{x + 2Tr) = —ip(x) and projecting out the states that are 
symmetric under ‘p{x) -G —ipix)^’^. A characteristic of 
the twisted sector of the orbifold boson is that it contains 
a number of primary helds whose scaling dimension is 
independent of Rq- Due to this construction the 0(2) 
symmetry of the free boson model is broken down to a 
Z 2 symmetry given by 'p(x) —>■ p>{x) + 'xRq, and the 
model retains its symmetry under ip{x) -G —ip{x) so that 
the natural symmetry group is Z 2 x Z 2 . 
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Figure 2. Each of the spin models is defined on a pair of 
parallel chains. The a operators act only on the red (top) 
chain and the r act only on the blue (bottom). The green 
region indicates a single site. Each of the models is then 
defined by a ID nearest-neighbour Hamiltonian. 


For a more detailed examination of these CFTs, in¬ 
cluding expressions for their spectra of primary fields, 
we refer to Ref. 5. 


B. Spin models at criticality 

In the following, we present three critical spin models 
whose thermodynamic limit is described by c = 1 CFTs. 
All three models will be defined on a line with periodic 
boundary conditions, and will possess an on-site Z 2 x Z 2 
symmetry. 

Each of these models have a line of critical points sep¬ 
arating an ordered phase from a disordered phase. For 
two of the models, the staggered XXZ model and the 
transverse field cluster model, the ordered phase displays 
symmetry protected topological order^®’^"*. In contrast, 
for the quantum Ashkin-Teller model the ordered phase 
displays conventional symmetry breaking of the Z 2 x Z 2 
symmetry. 

To clarify the on-site nature of this symmetry, we will 
assign two spin-1 /2 particles to each site, and colour them 
red and blue (Fig. 2). At site j, the Pauli spin operators 
for each spin are denoted cr“ and Tj for a,j3 = X, Y, Z. 


1. Staggered XXZ Model 

The staggered XXZ model is described by the Hamil¬ 
tonian 


N 


Hxxz = -J2 (^f 


i=i 

N-l 


-ZZ 




( 1 ) 




Note that this model has an 0(2) symmetry group, cor¬ 
responding to rotating each spin about the z-axis and a 
Z 2 symmetry corresponding to flipping each spin. 


In this work we will be mainly concerned with a Z 2 x Z 2 
subgroup of these symmetries that will remain on-site in 
all of our models, as we will incorporate this symmetry 
explicitly into our numerical MERA simulations. This 
symmetry is also important for considerations of symme¬ 
try protected topological order^. We can choose genera¬ 
tors for this subgroup as follows: 

N N 

-^1= n ^2 = n Y ■ (2) 

i=i i=i 

The generator corresponds to the spin-flip symme¬ 
try of the staggered-XXZ, and therefore to the ip ^ —(p 
symmetry of the free boson model. The generator S 2 
corresponds to rotating each spin by tt about the Z-axis 
and then flipping the spins. Thus it corresponds to the 
transformation (/?—)■ —ip + ir on the free boson. 

For /? = 0, A > — 1 the ground state corresponds to 
pairing the red and blue spins at each site in a maxi¬ 
mally entangled state, and thus has the structure of a 
product state. For /3 —>■ oo,A > —1 the ground state of 
this model pairs each blue spin with the red spin on the 
site to the right in a maximally entangled state. If the 
spin chain is defined on open boundary conditions there 
is a four-fold ground state degeneracy associated with 
unpaired spins at each end of the chain. This state is in 
the non-trivial symmetry protected topologically ordered 
state of Z2 X Z2 symmetry^. These two wave functions are 
the canonical representatives of the trivial and non-trivial 
(respectively) symmetry protected topologically ordered 
phases with this symmetry^®. 

At the phase transition f3 = 1 separating these two 
phases, this model is the well-studied XXZ model and is 
solvable via a Bethe Ansatz^®. This information can be 
used to identify how different values of A correspond to 
the parameter Rc that specifies the free boson model. 
The critical line defined by 

/? = !, A €[-72/2,1), (3) 

is known to be described by the free boson^’^ com- 
pactified on a circle of radius 

= -cos-^{-X). (4) 

TT 

Due to the exact solution, this antiferromagnetic model is 
used extensively in the study of quantum critical points 
and critical ground states^®^^°. We also note that the 
XXZ model with next-nearest neighbour interactions has 
recently been studied in the context of SPT phases^^. 
Moreover the study of phase transitions out of the SPTO 
phase of SO(S) symmetric systems in one dimension in 
Ref. 21 corresponds to the Heisenberg model in one di¬ 
mension and the value A = 1 in the free boson theory. 


2. Transverse field cluster model 

Cluster states are highly entangled many body states 
of spin-half particles that arose in the theory of quan- 
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turn computing. It is possible to simulate any quan¬ 
tum computation by making only local measurements 
on cluster states on an appropriate lattice^^. The clus¬ 
ter state corresponding to a linear arrangement of spins 
is the ground state of the following local Hamiltonian 
H = Recently, cluster state mod¬ 

els consisting of this Hamiltonian with various additional 
terms have been studied, from the perspective of quan¬ 
tum computing^^d9,20^ their phase structure as models 
with symmetry protected topological (SPT) phases^®’^^, 
natural models for the gapless edge physics of 2D SPT 
models^^’^^, and nonequilibrium dynamics arising in 
these models^®. 

We define the transverse field cluster model (TFCM) 
by the Hamiltonian 


N 

Htfcm = - y {<jf + rf + Aaf rf ) 




N-1 

j=i 


'ZXZ I zx z 
^3 ‘^3 + 1 ^3 ^ 3 + 1^3 + 1 




)■ (5) 


This model possesses an on-site Z 2 x Z 2 symmetry gen¬ 
erated by 


N N 

S, = l[af S2 = l[r^, ( 6 ) 

i=i i=i 

and possesses two distinct SPT phases for this symmetry. 
For /3 = 0, the ground state is a product state, i.e., in a 
trivial phase, whereas if /3 —>■ 00 the ground state is the 
well-studied cluster state, known to be within a Z 2 x Z 2 
symmetry protected topological phase^®. 

The transverse field cluster model can be obtained 
from the staggered-XXZ model through the transforma¬ 
tion 




^3 ^3 


^3 ^3 


(7) 






^3 ^3 ■ 


As this mapping is unitary, the ground state energy and 
phase transitions of the TFCM model are identical to 
those of the XXZ model. The mapping (7) preserves lo¬ 
cality, meaning that local operators confined to k sites 
of the staggered XXZ model map to local operators of 
the transverse field cluster model confined to the same 
k sites, and that this transformation can be performed 
using local unitaries applied to each site. Finally, this 
mapping respects the Z 2 x Z 2 symmetry of the models. 
As a result the ground states of the two models must be 
in the same phase^®. This is one way of seeing that the 
ordered phase of the transverse field cluster model pos¬ 
sesses symmetry protected topological order. Moreover 
we expect that the thermodynamic limit of TFCM on 
the line defined by Eq. (3) likewise corresponds to the 
boson with radius Rc determined by A as in Eq. (4). 


Finally the mapping (7) maps the symmetry operations 
Si , S 2 of the staggered XXZ model to the correspond¬ 
ing symmetry operators for the transverse field cluster 
model. Thus we can identify these two symmetry opera¬ 
tions in terms of the mappings tp —>■ —tp and tp —>■ —ip + n, 
respectively, in the free boson model. 


3. Quantum Ashkin-Teller model 


The quantum Ashkin-Teller (AT) modeR®’^^ is defined 
by the Hamiltonian 


N 

77 at = - E {^3 + ^3 + ^^3^3) 

Af-l 

i„xx I 


D ^3+1 x ‘j ‘j+l 


^xxx X 

^3 ^3 


( 8 ) 

)• 




This model consists of a pair of transverse-field Ising 
chains cr, r coupled by two- and four- spin terms. In par¬ 
ticular, when A = 0, we have a decoupled pair of Ising 
chains. 

This model also possesses an on-site Z 2 x Z 2 symmetry, 
generated by 

N N 

= S2 = X{rf (9) 

i=i i=i 


We note that with open boundary conditions this 
model is also unitarily related to the transverse field clus¬ 
ter model, but now through a nonlocal unitary transfor¬ 
mation 



( 10 ) 



Under this mapping the four-fold degenerate ground 
space of the transverse field cluster model maps to the 
four-fold degenerate groundspace of the Ashkin-Teller 
model. Although this mapping is nonlocal, in the 
bulk it maps local Z 2 x Z 2 -symmetric terms to local 
Z2 X Z 2 -symmetric terms, and thus can be viewed as a 
generalisation^^ of the map by Kennedy and Tasaki^®’^®, 
mapping the SPTO phase of the transverse field cluster 
model to the symmetry broken phase of the Ashkin-Teller 
model. Note, however, that boundary terms are trans¬ 
formed nontrivially by this map; local boundary terms 
can map to nonlocal ones. In particular, periodic bound¬ 
aries in the XXZ and TFCM models pick up extensive 
stringlike terms in the AT model. 

(Note that the Ashkin-Teller model also possesses a Z 2 
symmetry associated to swapping a ^ t, however this 
does not map to a local symmetry in the other models.) 

Because the models are unitarily related, the AT model 
will possess the same phase structure and ground state 
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Figure 3. Tensor network diagram for the scale invariant 
ternary MERA. Isometries are indicated by triangles (green), 
while disentanglers are rectangles (blue). The layered struc¬ 
ture accurately represents a wide range of length scales and 
captures the scale invariant physics of critical spin chains. 

energy as the transverse field cluster and staggered XXZ 
models (and, in particular, will possess the same line of 
criticality described by Eq. (3)). However, the phases 
that arise have quite different properties and the CFT 
corresponding to the critical line where j3 = 1 will have a 
different spectrum of primary fields. Based on finite size 
simulations^'^’"*^^ and CFT arguments^^’^^, the critical line 
of the AT model has been identified as the orbifold boson 
with radius 

Ro = Rc\ (11) 

where Rc is defined in Eq. (4). 

III. SCALE-INVARIANT MERA WITH 
SYMMETRIES 

Following the work of Refs. 44 and 45, we have inde¬ 
pendently developed code to optimise multi-scale entan¬ 
glement renormalisation ansatz (MERA)^^. This class 
of tensor networks can be used to efficiently describe 
ground states of critical quantum lattice models. Specif¬ 
ically, we give a brief introduction to the scale invariant 
ternary MERA, closely following Ref. 45, and describe 
the adiabatic crawling method which we used to stream¬ 
line the calculation of critical exponents at various lo¬ 
cations along the critical line of our c = 1 spin chain 
models. 


A. Elements of the MERA 

The ternary MERA shown in Fig. 3 is constructed from 
two kinds of tensors: isometries w 

Wi : Re+i —)■ wjwg = I^+i, (12) 


and unitary disentanglers u 

ui : Hf^ ^ nf, u\ui = If", (13) 

where T-in is the Hilbert space of dimension xn of one spin 
(one index) on layer i. 

Together these tensors perform a real space renormali¬ 
sation group (RG) transformation, with each layer of the 
structure corresponding to a description of the model on 
a different length scale. The free indices correspond to 
the degrees of freedom associated with the microscopic 
model of interest, and the isometries map 3 spins on 
layer £ onto an effective spin on layer £ -|- 1. The ap¬ 
proximation involved in this transformation is controlled 
by the bond dimension y; the dimension of the effective 
spin. If x^+i < xf, the full microscopic physics cannot 
be captured and a variational algorithm is used to en¬ 
sure the low energy subspace is retained. The maximum 
bond dimension used in the MERA is labelled x- The 
unitary tensors rearrange the local Hilbert space, locally 
removing entanglement and allowing x to be relatively 
smalb^. 

Generically, all tensors in the MERA may be differ¬ 
ent, however for translationally invariant states a sin¬ 
gle pair of tensors {ui,W£} characterise layer £. For 
scale invariant models, such as critical systems, the net¬ 
work becomes simpler still. The description of the model 
should be scale-invariant, and therefore all layers become 
identicab^ and the entire MERA is described by a single 
pair of tensors {u, w}, and all bond dimensions are x- 

The isometries and disentanglers are chosen to min¬ 
imise the energy of the spin model Hamiltonian using 
techniques that are described in Refs. 44 and 45. This 
results in an iterative algorithm where each step updates 
either the isometries or the disentanglers and each step 
can be performed at a cost 0(x®), although this can 
be reduced by modifying the network to include further 
approximations, such as spatial symmetries^® or on-site 
symmetries^®. 


B. C/-symmetric tensors 

The incorporation of symmetries can decrease the re¬ 
sources required to optimise the MERA. If a model pos¬ 
sesses an on-site symmetry, such as the Z 2 x Z 2 symmetry 
occurring in the models described in section H, these in¬ 
ternal symmetries may be exploited to further reduce the 
resource requirement by using Q-symmetric tensors'^^. 

A tensor T^^aV ,^-symmetric if it is 
invariant under the action of the group Q on each index. 
For example, the requirement for the isometry described 
above to be ^-symmetric is as follows 
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Ua V, 


where U, V, W, X are unitary representations of Q. 

By Schur’s lemma, the tensor decomposes into blocks, 
where each block transforms as one of the irreducible rep¬ 
resentations (irreps) of In general, due to the block 
structure of the tensors, we can decompose the indices 
a —> {c,d), where c labels the irrep and we will call it 
the charge index, and d is the degeneracy index^®. The 
charge structure is completely specified by the group G- 
The degrees of freedom for the particular model of inter¬ 
est are completely described by the index d. 

In our case we have the group G — Z 2 X Z 2 = {x, y\x^ = 
= e,xy = yx} whose irreps are given by the character 
table 



e 

X 

y 

xy 

-0(1.1) 

1 

1 

1 

1 

D(i.-i) 

1 

1 

-1 

-1 

D{-is) 

1 

-1 

1 

-1 

O’(-i.-i) 

1 

-1 

-1 

1 


where the Dc are the irreps and the possible charges 
c = ( 51 , 92 ) are (1,1), (1,-I), (-1,1), (-1,-1). Since 
the group is Abelian, tensor products of irreps results 
in a new irrep that can be obtained by elementwise 
multiplication in the character table. For example, 
Z> 13(1,- 1 ) = This product of taking 

tensor products results in multiplicative operation on 
charges ci (g) C 2 = C 3 that is often called the fusion rule. 
Using this fusion rule, we can construct higher order ten¬ 
sors from vectors. Condition (14) becomes a charge con¬ 
servation rule"^^. In numerics it is possible to identify 
the nonzero blocks of a given tensor by checking that 
the block conserves charge. The condition for a nonzero 
block Bi is 


Cn+1 Cn-\-m 

\ ••• I 

B. ^ X{c,= Cl. (15) 

£GOut 

Cl C2 Cfi 

Using these block tensors provides us with three main 
advantages in optimising the network and extracting 
physical information. 



1. Reduction in variational parameters 

A large number of the blocks are fixed to be zero by the 
symmetry, so do not have to be stored or manipulated 
during optimisation. This decrease in variational param¬ 
eters leads to a decrease in the storage space needed for 
the network. It also leads to a multiplicative decrease 
in the computational time (ay® —>■ in the case of 

Z 2 X Z 2 ) allowing either larger y to be accessed or a 
decrease in the overall runtime. 


2. Selection of symmetry sector 

By decomposing the tensors into blocks, any operator 
applied to the spin chain can be classified according to 
how it transforms under the symmetry. A operator Oq^ ,^3 
has charge ( 51 , 52 ) under the Z 2 x Z 2 symmetry if 

SiOS\ = 5 iO, S 2 OSI = 52 O, (16) 

where ^i, S 2 are the symmetry generators described in 
section II. Operators carry a unique charge, so operators 
of the form O = aiOi i + a20i._i are forbidden. 


3. Nonlocal operators 

As we know the action of the symmetry on the tensors 
{u,w}, we can trivially compute expectation values for 
a class of highly nonlocal operators, with computational 
cost which barely exceeds that of local operators^®. This 
is discussed below with respect to scale invariant opera¬ 
tors and conformal data. 


C. Optimising a MERA for a spin model 

Given a local, critical spin Hamiltonian H = hj , a 
MERA description of the ground state \'if) can be gen¬ 
erated by optimising the tensors. Here we briefly review 
the optimisation algorithm, following Ref. 45. 

The MERA for a critical model is built from two sec¬ 
tions. First, a number (usually 2-3) of transitional layers 
are used, which are translationally invariant but not scale 
invariant, with tensors {ue,we}. These allow the bond 
dimension of the remaining (scale invariant) tensors to 
be chosen independently of the physical dimension of the 
spins. Generically, the Hamiltonian H will contain RG 
irrelevant terms, breaking scale invariance. Since these 
terms become suppressed at larger scales, the transitional 
layers reduce their effect. 

Above the transitional layers, the MERA is built from 
a unique pair of tensors {m,?c}, characterising the scale 
invariant nature of the critical model. We allow the di¬ 
mension of the upper and lower indices of u to differ; it 
is not unitary, but Eq. (13) holds. This preserves the es¬ 
sential structure of the MERA, but allows for increased 
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Figure 4. The descending superoperator. This computes the 
reduced density matrix of the model at level t of the RG. Just 
one of the contractions which must be performed to optimise 
the MERA tensors. The cost is 0{x^), however this can be 
reduced as discussed in the text. 


values (scaling dimensions for relevant operators and cen¬ 
tral charge) is defined as v = Cj — Cj-i, where Cj is 
the conformal data on iteration j. When this falls be¬ 
low some threshold, the point is declared converged, and 
the algorithm is repeated on the next point, using the 
previously converged tensors as a start point. If the vari¬ 
ation is sufficiently small, the ground state of the new 
Hamiltonian is very similar to the previous ground state, 
and few convergence iterations are required. Typically, 
convergence of the remaining 60-100 points can be com¬ 
pleted in a time similar to that required to converge the 
initial point. 


numerical efficiency^^. A pair of bond dimensions, xi 
Xu on the lower and upper indices of u characterise the 
MERA. 

To initialise the MERA, a pair of bond dimensions 
iXhXu} is chosen, and random tensors generated. A sin¬ 
gle transitional layer is initialised to allow xi to differ 
from the physical dimension d. These tensors are then 
iteratively optimised to minimise the energy {tjj\ H \ tp) = 
Tt{H |'!/')('0I) as discussed in Refs. 44 and 45. 

Optimising the MERA requires many contractions of 
networks analogous to those shown in Fig. 4, where is 
the reduced density matrix describing the system on layer 
£. This contraction can be computed in time O(xfx^). In 
numerical implementation, Xu < Xi allows for decreased 
runtime without apparent degradation of the results. By 
assuming that most of the eigenvalues of the reduced den¬ 
sity matrix p are small, we can modify the network to im¬ 
prove the scaling as follows"*®. We assume we can choose 
a projector P — vv^ such that p = Pp -I- e, where e is 
a small error. Then x is the rank of the projector. In 
the networks, p is replaced with Pp, giving the improved 
efficiency at the expense of optimising the tensor v. Nu¬ 
merical evidence shows that x may be chosen 0{x) rather 
than 0{x^), thus providing a scaling of 0(x^x). 

As the optimisation proceeds, the change in the ten¬ 
sors between iterations decreases. Once the network is 
changing sufficiently slowly, a new transitional layer is 
added by promoting the lowest level of the scale invari¬ 
ant portion. This process is repeated until it does not 
produce significant improvement of the ansatz (usually 
once there are 2-3 transitional layers). 


D. Adiabatic crawling 


E. Extracting conformal data from the MERA 

The properly converged MERA provides a compact 
approximation to the ground state, and we can extract 
physical properties of the GET that describes the spin 
lattice model at criticality by optimizing a MERA. Here, 
we describe how to obtain the conformal parameters from 
the tensors in the MERA, following Refs. 12 and 13. 

One part of the conformal data for a GET is the spec¬ 
trum of scaling dimensions associated with a primary 
field 4>- These are the eigenvalues of the rescaling opera¬ 
tor in the field theory. In the MERA, the isometric tensor 
w performs a rescaling operation on the spin model, so 
allows us to extract the scaling dimensions*^. By finding 
the operators which are fixed under the one site scaling 
superoperator 



(17) 


we find the scaling operators and their dimensions. The 
eigenvalues of S are then related to the scaling dimen¬ 
sions via = log 3 (A 0 ). 

We can also compute the scaling dimensions for a class 
of nonlocal scaling operators by making use of the sym¬ 
metry. Nonlocal scaling operators take the form 


We are simulating families of models with properties 
that vary continuously along a critical line, and so we 
can make use of the solution from a converged MERA 
at one point as a starting point to speed up convergence 
at a neighbouring point. We call this approach the adia¬ 
batic crawling technique. First, the MERA is converged 
from random tensors at a start point somewhere along 
the line. The velocity rate of variation of the conformal 


N-l 

0 = {Ugh<»ON, (18) 

j=-cx> 

where Ug is a unitary representation of Q. Since the sym¬ 
metry operators commute with the MERA tensors, ap¬ 
plying one layer to an operator of this form and utilising 
Eq. (12-13) gives a new operator of the form Eq. (18). 








Eigenoperators of the g-twisted scaling superoperator 


A 



allow extraction of the conformal spectrum of the CFT 
arising from taking the thermodynamic limit of the spin 
model with boundary conditions which are twisted by 
the group element g; see Ref. 13. Setting g — e, S is 
recovered, giving the spectrum on periodic boundaries. 

The spectrum of Sg can be extracted from the con¬ 
verged MERA for negligible additional cost simply by 
diagonalizing 

The central charge of the theory can be extracted from 
the fixed-point density matrix using the scaling of the 
entanglement entropy^^ 

Q 

‘^'critical — ^ /^) ~b ^5 

where c is the central charge of the CFT, L is the block 
length, a is the lattice scaling (here, a = 1 by definition) 
and k is some constant. The one-site reduced density 
matrix is obtained by symmetrising over the two ways 
of tracing out one site of the fixed-point density matrix 
P 2 ^ (Fig. 4). The central charge is then obtained by 

c = 3 (s' - S j ) , (21) 

where S is the von Neumann entropy of the density ma¬ 
trix. 


IV. NUMERICAL RESULTS 

In this section, we present results obtained from a 
Z 2 X Z 2 symmetric MERA with {xhXu} = {20,12} and 
X = 60 (see Sec. IIIC) applied to the three models de¬ 
scribed in Sec. II. Three transitional layers are used to 
ensure scale invariance of the Hamiltonian. We use the 
adiabatic crawling technique, starting at = 1.25 and 
using a velocity threshold of = 9 x 10“®. Both ground 
state energies and conformal data are presented along the 
critical lines of the three models given by Eq. (3). The 
conformal spectrum is shown to be consistent with the 
identifications made in Sec. II. We also compare these 
MERA simulations with the result of exact diagonalisa- 
tion studies for the models of interest. 


A. Ground state energy 

As a result of the Bethe Ansatz solution for the XXZ 
model, the ground state energy is known exactly for all 


-2 


-2.5 


O 

-3 


-3.5 


X- 

X 

X AT 

+ TFCM 

XXZ - 

- Exact 

X 

%■ 

% 

% 

_1_1 

X. 


0.5 1 1.5 2 


CFT Radius 



Figure 5. Ground state energies (GSE) extracted from the 
MERA and their relative errors AE. For the bond dimension 
used here, these remain 0{1Q~*) as the Hamiltonian parame¬ 
ter A is varied. The numerical GSE remains an upper bound 
on the true GSE. Note that the CFT radius Rc is related to 
the spin model coupling parameter via Eq. (4) 


three models considered here, and this provides a useful 
benchmark; see Fig. 5. The relative error in the ground 
state energy is given by AE = (Aexact - AMERA)/Aexact- 
The exact solution is obtained by numerically integrat¬ 
ing the Bethe Ansatz solution^®. The MERA is an ex¬ 
plicit wavefunction, and as such, its energy cannot be less 
than the true ground state energy; thus AE > 0. The 
ground state energy obtained using the MERA for all 
three models is consistent with the exact solution, with 
relative error positive and of order 10“^ for the full range 
of A considered; see Fig. 5. 


B. Conformal data 


We now present the conformal data obtained from the 
MERA for our three models. We show that this data 
is consistent with the CFTs identified as the thermody¬ 
namic description, namely the Z 2 -orbifold boson CFT 
for the Ashkin-Teller model; and the boson CFT as¬ 
sociated with both the XXZ and transverse field cluster 
models. 
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Figure 6. Central charges extracted from MERA simulations. 
The CFTs associated with these models are expected to be 
c = 1 (shown as solid line) for the full range of Rc considered 
here. 



Figure 7. The lowest scaling dimensions (A^) of the (1, —1) 
charge sector of the three models. The XXZ and transverse 
field cluster models agree well with the boson, whilst the 
AT model is consistent with the fixed scaling operator from 
the orbifold theory. Note that the and orbifold CFT radii 
are related via Eq. (11). 
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1. Central Charge 

One of the pieces of data required to specify a CFT 
is the central charge. It labels classes of CFT and is 
identically c = 1 for all values of R in both the boson 
and orbifold boson CFTs. The values obtained from the 
MERA simulations are shown in Fig. 6. These are usually 
within 2% of the expected value c = 1 for the full range 
of Rq. We note an increased deviation from c = 1 at the 
left (Rq ^ 0.5) end of the critical line. Here, helds in 
the Z 2 X Z 2 symmetric sector are crossing = 2, the 
RG relevant/irrelevant threshold. Some of these helds 
are not symmetric under the 0(2) symmetry discussed in 
section iP®. Near R^ = 1.6, we also see a large deviation 
from c = 1 for the TFCM data. This region occurs close 
to a crossing of two helds in the Z 2 x Z 2 symmetric sector 
of the theory. Again, this crossing does not occur in 
the 0(2) symmetric theory^®. 

Our simulations approach, but do not include, the end 
point A = 1, {Rq = 2). At this point, the XXZ model 
becomes the spin-1/2 Heisenberg model and the line of 
criticality terminates. As we approach this point, the 
CFTs have an irrelevant primary operator with a scal¬ 
ing dimension that is decreasing as Rc increases and be¬ 
comes marginal precisely at Rq = 2. This operator leads 
to corrections to the scale invariant behaviour that, as 
has previously been noted"*®, complicate studies of this 
model using MERA simulations. 

This limiting case could be further investigated in two 
ways. One approach would be adding more transitional 
layers to the MERA, which is expected to move the 
Hamiltonian towards the RG fixed point. However, many 
layers may be required because this term only decays 
logarithmically. Alternatively, modifying the model by 
adding hnely tuned terms (for example, a next nearest 
neighbour coupling) could remove the marginally irrele¬ 
vant operator without changing the continuum limit^®’^®. 
Both of these approaches are beyond the scope of the cur¬ 
rent investigation. 


2. Scaling dimensions 

We now turn to the scaling dimensions extracted 
from the MERA. These scaling dimensions can be classi¬ 
fied according to their Z 2 x Z 2 symmetry sector. As the 
parameters of the models are varied along the lines of crit¬ 
icality, the scaling dimensions of the two CFTs can vary 
continuously. As noted in section H, the orbifold CFT 
contains a fixed sector of primary fields whose scaling di¬ 
mensions do not change with Rq , whereas for the free bo¬ 
son CFT all scaling dimensions (other than the identity 
and related operators) will vary. The behaviour of these 
scaling dimensions provides a good pointer as to the type 
of c = 1 CFT that corresponds to each of the spin models. 
In Fig. 7, we focus on fields with small scaling dimension 
in the (1,-1) symmetry sector. Previous evidence*®’*® 
indicates that the accuracy of the MERA decreases with 
increasing scaling dimension, so it is expected that the 
agreement will be better in this region. The distinction 
between the CFTs is already evident, and the scaling di¬ 
mensions obtained from the XXZ and TFCM models are 
consistent with the S'* theory, and the AT model has the 
expected fixed scaling dimension present in the orbifold 
theory. 

We next examine the full relevant spectrum of the scal¬ 
ing operators for all three models, up to a maximum value 
of A^ = 2. All sectors of the XXZ and TFCM models are 
expected to have the same conformal dimensions; those 
of the S* theory. Close agreement is seen in the MERA 
simulations of the two models (Fig. 8 a,b), although the 
accuracy deteriorates for larger scaling dimensions. The 
conformal held theory has an inhnite number of scal¬ 
ing operators, however the one-site scaling superopera¬ 
tor used to obtain has a hnite number of eigenvalues. 
This results in only a hnite subset of the dimensions be¬ 
ing recoverable, with the smallest being most accurately 
reproduced. In particular, A^ = 0 in the (1,1) sector 
(Fig. 9) corresponds to the identity operator. This seal- 
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Figure 8. The behaviour expected from the 5"^ CFT is recovered by the XXZ and TFCM (a,b), whilst AT simulations recover 
the orbifold behaviour (c,d). The accuracy of the recovered scaling dimensions degrades for larger A^. Note that the 
and orbifold CFT radii are related via Eq. (11). 


ing dimension is recovered perfectly due to the isometric 
constraint on the tensor w. 

The AT model is expected to reproduce the orbifolded 
theory, with the signature fixed sector being equally split 
over the (1, —1) and (—1,1) charge sectors. We see this 
indeed occurs in Fig. 8c, with good agreement for the 
lower scaling dimensions. 

In the (1,1) sectors (Fig. 9), we see the deviation in¬ 
crease as we move towards = l/i?o = 2. This also 
occurs in the XXZ and TFCM models. Recall that, at 
this endpoint, all three models are unitarily equivalent to 
the Heisenberg model. 


C. Exact diagonalisation of the quantum 
Ashkin-Teller model 

It is well known that the low energy spectrum of finite 
size spin chains should correspond to the spectrum of 
primary and descendant fields of the associated CFT®. 
As such the results of exact diagonalisation (ED) with 
periodic boundary conditions can be directly compared 
with those of both CFT and the MERA simulations. Eor 
a review of the ED technique, see for example Ref. 50. 

As a further test of the MERA simulations, we per¬ 
formed ED calculations for the Ashkin-Teller model with 
periodic boundary conditions for chains up to L = 12 in 


length. At this length, finite-size effects are still present 
but we are nevertheless close to full convergence. As can 
be seen in Fig. 9a, the gap is closing like 1/L, as is ex¬ 
pected of a critical model®. 

The spectrum of the quantum spin chain can be related 
to that of the CFT once the ground state energy (Eg) and 
the overall normalisation of the Hamiltonian are chosen 
appropriately®^. In our problem we need to choose a 
normalisation of H for each and we achieve this simply 
by rescaling the gap {Ei — Eq) of the spin model. Then 
the rest of the conformal spectrum can be inferred from 
the spectrum of the spin chain as follows 

Ak = LxaiX)iEk-Eo), ( 22 ) 

where, anticipating the orbifold CFT, the a{X) is cho¬ 
sen to fix the gap Ai = 1/8. This value corresponds to 
the scaling dimension of the operators with fixed scal¬ 
ing dimension in the (±,t) sectors (Fig. 8d). Thus the 
first two scaling dimensions of the CFT are correct as 
a result of our choice of normalisation and the zero of 
energy. Obtaining the remaining levels represents a con¬ 
firmation of the CFT identihcation. The same approach 
was taken in previous studies of the quantum Ashkin- 
Teller model"*^®’^^ , although there an analytic form for a 
was proposed. 

We see good qualitative agreement between the two 
methods. A selected Z 2 x Z 2 symmetry sector is pre- 
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Figure 9. Results of an exact diagonalisation study of the Ashkin-Teller model in the symmetric (1,1) sector, a) The energy 
gap for a range of values of Ro and system sizes L. The gap closes as we approach the thermodynamic limit. The gap is 
plotted for various I/Rq values; in blue (x) is the point 1/Ro ~ in magenta (+) the point l/i?o ~ 2. Lines represent a 

linear fit to the ED data, b) A comparison of the finite-size scaling results for L = 12 with the MERA and the GET, for the 
Z 2 X Z 2 symmetric sector. 


sented in Fig. 9 b; we see comparable agreement for all 
other sectors. The agreement between the two methods is 
extremely good for lower scaling dimensions but becomes 
less so for higher ones above « 1.5. The ED and CFT 
data match very precisely at the point l/i?o = 1 for all 
scaling dimensions up to ~ 2 , where the model is 
equivalent to two uncoupled critical Ising models. Hence, 
away from the decoupling point we interpret the decrease 
in precision with increasing scaling dimension as a finite- 
size effect rather than a result of the Lanczos-based ED 
algorithm. 

It is interesting to note that close to the point 1 /= 
2 the MERA and the ED results agree with each other 
rather better than either one agrees with the CFT values. 
This is the limit in which the finite size corrections to the 
eigenvalues for the spin chain, which are well understood 
for ED52, increase in size due to an irrelevant primary 
field whose scaling dimension is approaching 2. It would 
be an interesting avenue for further study to understand 
whether the errors in scaling dimensions obtained from 
MERA simulations can be understood in a similar way 
to those of ED. 

We note that finite-size scaling has previously been 
performed for the Ashkin-Teller chain, in order to iden¬ 
tify the primary fields for periodic and anti-periodic 
boundary conditions^^’"^^’®^. To the best of our knowl¬ 
edge, our ED calculations use a system size equal to the 
largest used in all previous studies^°. Here we have ex¬ 
tended those calculations over a larger range of 1 /as 
well as extracted scaling dimensions for all RG relevant 
fields. 


D. Nonlocal operators and twisted boundaries 

Another class of scaling operators present in the sym¬ 
metric MERA are the nonlocal operators involving a half 


infinite string of symmetry operators terminating in a lo¬ 
cal operator. The scaling dimensions of these operators 
in the CFT arise in ED studies of the corresponding spin 
model with boundary conditions twisted by a group el¬ 
ement. Recall that the scaling dimension of these non¬ 
local operators can be obtained from MERA simulations 
at negligible additional computational cost. 

In Fig. 10, we show these nonlocal scaling operators in 
the ( 1 , 1 ) symmetry sector for two kind of twist (corre¬ 
sponding to the two of the group elements of Z 2 x Z 2 ) 
obtained from the MERA simulation of the Z 2 x Z 2 sym¬ 
metric AT model. Note that Fig. 10 and Fig. 8 show 
different data. The orbifold CFT is such that the lo¬ 
cal (untwisted) spectrum in the four symmetry sectors is 
equivalent to the twisted spectrum in the ( 1 , 1 ) sector. 
This equivalence is apparent by comparing Fig. 8 c (d) 
with Fig. 10b (a). The deviation of the scaling dimen¬ 
sions of these nonlocal fields from the CFT expectation 
appears to be comparable to that of the local fields. 


E. Observations and Comments 

All the models considered here have a larger on-site 
symmetry than the enforced Z 2 y- Z 2 . In particular, the 
AT model is invariant under swapping a t, leading 
to a nonabelian ZI 4 invariance. The symmetry group of 
the staggered XXZ model contains an SO{2) symmetry, 
for SO{2) rotations of all spins about the z-axis, and a 
Z 2 symmetry corresponding to 7r-flips around the x-axis, 
just as the free boson CFT as described in Sec. H. This 
maps to an on-site symmetry in TFCM, but becomes a 
nonlocal symmetry in the Ashkin-Teller model unless the 
50(2) angle is tt (in which case this symmetry reduces 
to the Z 2 X Z 2 symmetry we have enforced). 

Enforcing only a proper subgroup of the full symmetry 
means that the results of the MERA simulations do not 
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Figure 10. Scaling dimensions of nonlocal operators extracted 
from MERA simulations. These correspond to the spectrum 
of spin chains with boundary conditions twisted by the group 
elements defined in Eq. (9). 


respect the full symmetry of the model. We note that 
this leads to avoided crossings in the conformal dimen¬ 
sions extracted from the simulations, as seen in Fig. 11. 
The two operators in this plot transform like different 
representations of and so are forbidden to mix. Un¬ 
der the Z 2 X Z 2 subgroup used in our simulations, such a 
mixing is allowed, and is indeed observed. We expect this 
gap to close as the approximation is improved, leading to 
an approximately enforced ZI 4 . 

We also note that under the full conformal symmetry, 
all primary fields transform like different irreps, meaning 
they are all uncoupled. As the MERA does not incor¬ 
porate the full (continuous) conformal symmetry in its 
structure, any MERA-based numerical method cannot 
keep all such fields from mixing, and as such we expect 
avoided crossings to be observed even if the full on-site 
group is enforced. 

We also note that the choice of starting point for our 
adiabatic crawling method appears to have an effect on 
the accuracy of the resulting converged MERA. Specifi¬ 
cally, at values of Rc or Rq where previously irrelevant 
operators cross over to become relevant, increased errors 
arise that slowly die away as the simulation proceeds. A 
similar situation occurs when the scaling dimensions of 
two relevant operators cross as Rc and Rq are varied. 
Crawling from multiple points, and stopping at crossings 



Eigure 11. Only the Z 2 x Z 2 subgroup of the U 4 symme¬ 
try of the Ashkin-Teller model is enforced in the MERA. 
This can lead to coupling between scaling operators which 
one would expect to be forbidden, in turn leading to avoided 
crossings in the conformal dimensions. This data has xi = 
20 = x/4, = 16 and 4 transitional layers. 


and when new operators become relevant may reduce this 
behaviour. 


1. Choice of X 

In this work, we have made a choice of the bond di¬ 
mensions {xiiXu\ = { 20 , 12 } and % = 60 but we have 
investigated a range of x values both higher and lower 
than this. These values were set by our available com¬ 
putational resources but we do not expect the essential 
conclusions to be altered by larger values of y. The scal¬ 
ing of the error in the MERA energy was investigated in 
Ref. 45 and elsewhere and is one of the few quantitative 
methods of studying convergence of MERA simulations 
reported in the literature. As discussed above, we can 
compare the energy of our MERA simulation with the 
exact result available from the Bethe ansatz solution and 
we also observe that for increased x the error in the en¬ 
ergy is reduced, as one would expect. Indeed the error in 
Fig. 5 is essentially the same size as reported in Ref. 45 for 
the staggered Heisenberg and XX models and the same 
value of X so we believe that our simulations have compa¬ 
rable accuracy to other implementations of MERA with 
comparable values of x- 

We have also qualitatively studied the convergence of 
the smaller scaling dimensions and central charge to the 
values predicted by CFT. The error in both decreases 
slowly with increasing x, with the smallest scaling dimen¬ 
sions being recovered more accurately than those higher. 

The qualitative features we see in the spectra of scal¬ 
ing dimensions are also robust and observed over a range 
of choices of x- For example the avoided crossings 
remain present even when increasing to {xijXu,x} = 
{20, 20,80}. Finally, we note that the deviation between 
the results obtained from the MERA and the CFT result 
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as we approach the endpoint of the critical line corre¬ 
sponding to = 1/ = 2 continues to persist even as 

X is increased. 

V. SUMMARY AND DISCUSSION 

In this work, we have used independently-developed 
scale invariant Z 2 x Z 2 symmetric MERA code to in¬ 
vestigate critical quantum spin lattice models with this 
symmetry. Specifically, we have used this to simulate 
the staggered XXZ and transverse held cluster models, 
extracting conformal data consistent with the boson 
conformal held theory being the thermodynamic descrip¬ 
tion of both spin chains. This c = 1 theory has a pa¬ 
rameter Rc which can be varied, leading to continuously 
varying critical exponents, a behaviour which has been 
replicated in the MERA. In addition, we have extracted 
conformal data for the Ashkin-Teller spin chain, identify¬ 
ing the .Z 2 -orbifold boson as the appropriate GET. Once 
again, the behaviour of the scaling helds with the param¬ 
eter Ro has been recovered, including the fixed sectors 
which are signatures of this theory. 

We have introduced a crawling method that allows for 
efficient optimisation of MERA along the critical lines. 
By examining the symmetries of the three models, we 
have identified some of the limitations of this technique, 
and how enforcing only a subgroup of the full symme¬ 
try is revealed in the conformal dimensions. We have 
also identified similarities between the numerical results 


of MERA and exact diagonalisation, in the behaviour of 
errors at different regions along the critical line. 

We note that the staggered XXZ and transverse field 
cluster models possess phases with nontrivial Z 2 x Z 2 
symmetry protected topological order. Such phases sup¬ 
port gapless edge modes that are protected by the sym¬ 
metry, a property not shared by SPT trivial phases^. Un¬ 
derstanding the critical theory occurring at the transition 
between the trivial and SPTO phases may provide insight 
into the fate of these edge states at the phase transition 
and their properties within the phase^^. We have seen 
how tuning parameters in the spin models leads to con¬ 
tinuous variation in the critical theory, despite the fact 
that we are investigating the transition from a single SPT 
phase to a single trivial phase. 

Recent developments of MERA allow for the incor¬ 
poration of conformal defects, including interfaces and 
boundaries®"^’^®. These numerical tools offer the possi¬ 
bility to study the gapless edge modes via the ‘domain 
wall boundary conditions’ described in Ref. 21 as well as 
interfaces between different SPT phases^®. 
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